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ABSTRACT 

We present results from A^-body + magnetohydrodynamical simulations of 

CIh- merging clusters of galaxies. We find that cluster mergers cause various char- 

O ■ acteristic magnetic field structures because of the strong bulk fiows in the intr- 

^ ■ acluster medium. The moving substructures result in cool regions surrounded 

, ^, - by the magnetic field. These will be recognized as magnetized cold fronts in 

the observational point of view. A relatively ordered magnetic field structure is 

j> ■ generated just behind the moving substructure. Eddy-like field configurations 

^ ■ are also formed by Kelvin-Helmholtz instabilities. These features are similarly 

t^ ■ seen even in off-center mergers though the detailed structures change slightly. 

. . The above-mentioned characteristic magnetic field structures are partly recog- 

Q . nized in Faraday rotation measure maps. The higher absolute values of the 

OO . rotation measure are expected when observed along the collision axis, because of 

the elongated density distribution and relatively ordered field structure along the 



axis. The rotation measure maps on the cosmic microwave background radiation, 
which covers clusters entirely, could be useful probes of not only the magnetic 
field structures but also the internal dynamics of the intracluster medium. 

Subject headings: galaxies: clusters: general — magnetic fields — magnetohy- 
drodynamics: MHD — X-rays: galaxies: clusters 



1. Introduction 

Clusters of galaxies have plenty of hot plasma as well as galaxies and dark matter 
(DM), which is called intracluster medium (ICM). Some observational evidence indicates that 
ICM is magnetized. For example, some clusters have diffuse non-thermal synchrotron radio 



emission that is called radio halos or relics, whic h shows that there are magnetic fields as well 
as rel ativistic electrons in the intracluster space (JGiovannini et al.lll999l : iKempner fc Sarazin 
2OOII ). In addition, Faraday rotation measure observations of polarized radio sources such as 
radio lobes and A GNs behind and /or in clusters indicate that the mag:netic fie ld structures 
are quite random (IClarke et al.ll200ll : IVogt fc EnsslinI l2003l : iGovoni et al.ll2006l ). 



Comparing the synchrotron radio flux with the hard X-ray one (or its upper limit) 
due to the inverse Compton scatterings of cosmic microwave background (CMB) photons, 
we are able to estimate the volume averaged magnetic field st rength (or its lower limit) . 



Typi cally, the strength of ^ 0.1/xG is obtained with this method (JFusco-Femiano et al. 



1999. 



2005h though those detections of non-thermal hard X -ray components are still controversial 
(JRossetti. fc Molendill2004l : lFusco-Femiano et al.ll2007l ). On the other hand, somewhat higher 
values (^ a severa l /iG) tend to be obtained with the Fara day rotation measure method 



(IClarke et all I2OOII : IVogt fc EnsshnI l2003l : iGovoni et all I2OO6I ) though these results depend 
on the detailed mag netic field structures that are not fully understood (JEnsslin fc Vogtll2003l : 
Murgia et alll2Q04l ). 



Although the magnetic energy density is typically ^ a few percents of the thermal one 
in the intracluster space, it is believed that the magnetic fields play a crucial role in various 
aspects of ICM. It is likely that non-thermal particles are accelerated via shocks (ISarazin 



1999 



Takizawa fc NaitQll200d: IXotani fc Kitayama 



2000 



Miniati et al 



2001 



Rvu et al 



2003 
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2003; 



2OO1 
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Inoue et al. 



Qhno et al. 



2005h and /or turbulence (|Rolandl 11981 
2OO2I : iFujita et al.lbood : ISrunetti et al. 



1987: 



200 



I. iSchlickeiser et al. 

04). where the magnetic 



fields play a definitive role in most theoretical models of the particle acceleration. Heat con- 
duction is probably depends on the magnetic field configurations because charged particles 
cannot freely move in the direction perpendicular to the field lines. As we wrote before, 
magnetic field seems to be a minor component in global ICM dynamics. However, it is pos- 
sibly important in relatively small scales, where fluid instabilities such as Ray leigh- Taylor 
and Kelvin-Helmholtz ones might be suppressed by magnetic tension. 

Cluster mergers have a significant impact on ICM magnetic field evolution. Turbu- 
lent motion excited by mergers would amplify the field strength via dynamo mechanism. 
Moving substructures are e xpected to sweep the fi e ld lines and form the c old subclumps sur - 



rounded by the field lines (jVikhlinin et al. 



2OOII). lAsai et all (120041 ) and lAsai et all (120071 ) 



performed two and three dimensional magnetohydrodynamical(MHD) simulations of moving 
cold subclumps in hot ICM in rather idealized situations, respectively, which confirmed that 
expectation. This field structure might be responsible for the suppression of heat conduction 
and fiuid instabilities, which is essential for the maintenance of cold fronts. 



A lot of numerical simulations about merging clusters have been done, most of which 



are A^-body + hydrodynamical simulations (IRoettiger et all Il996l: iTakizawal Il999 



Ricker fc Sarazin 



200 



20061 : iMcCarthv et al. 



Ritchie fc Thomas 



2002; 



Ascasibar fc Markevitch 



2006; 



2000 



20071 : ISpringel fc Farrarl 120071 ) . Although these simulations gi 
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ive us 

a great deal of understanding of the structures, evolution, and observational implications 
for merging clusters of galaxies, they make only a limited contribution to investigate the 
magnetic field structures. MHD simulations are essential in this r egard. However, A^-body 
+ MHD simulations are rather rare though IRoettiger et al.l (119991 ) did pioneering work. It 
is true that LaRr auRian particle methods based on smoothed particle h ydrodynamics (see 



1999. 



Monaghan 1 1 19921 ) are extensively used in cosmological MHD simulations (iDolag et al. 
2OO2I ). Considering that Eulerian mesh codes are essentially better at following evolution 



of magnetic fields, simulations based on such codes are highly desirable. In this paper, we 
present the results from A^-body + MHD simulations of merging clusters of galaxies and 
investigate characteristic magnetic field structures during mergers and their implications. 



The rest of this paper is organized as follows. In % 
methods and initial conditions for our simulations. In 
summarize the results and discuss their implications. 



I we describe the adopted numerical 
13] we present the results. In ^ we 



2. The Simulations 



2.1. Numerical Methods 



Numerical methods used here are basically similar to those in ITakizawal (120061 ) except 
that ideal MHD equations are calculated for ICM instead of ordinary hydrodynamical ones. 
In the present study, we consider clusters of galaxies consisting of two components: colli- 
sionless particles corresponding to the galaxies and DM, and magnetized gas corresponding 
to the ICM. When calculating gravity, both components are considered, although the for- 
mer dominates over the latter. Radiative cooling and heat conduction are not included. In 
calculating the M HD equations for t he ICM, we use a simple linearized Riemann solver origi- 
nally proposed bylBrio fc Wul (119881 ) and introduced into astrophysical MHD simulations by 
Rvu fc JonesI (119951 ). This Riemann solver is a similar and simpler version of Roe^s me thod 
( JRoelll98ll ). Using the MUSCLE approach and a minmod TVD limiter (see lHirschlll990l ). we 
obtain second-order accuracy without any numerical oscillations around discontinuities. To 
avoid negative pressure, we solve the equations for the total energy and entropy conservation 
simultaneously. This method is oft en used in astrophysical hydrodynamic si mulations where 
high Mach number flow can occur (JRyu et al.l 1 19931 : IWada fc Normanll200ll ). 



Gravitational forces are calculated by the Particle-Mesh (PM) method with the standard 



Fast Fourier Transform technique for the isolated boundary conditions (see lHockney fc Eastwood 
19881 ). The size of the simulation box is (9.44 Mpc)^. For typical runs, the number of the grid 
points is (256)^, and total number of the A^-body particles is (128)^, which is approximately 
2.1 X 10^ 



2.2. An Equilibrium Cluster Model 



W e consider mergers of two virialized subclusters with an NFW density profile (JNavarro et al. 
1993) in the ACDM universe (^o = 0.25, Aq = 0.75) for DM, 
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where 6c ^ Pcrit, and r^ are the characteristic (dimensionless) density, critical density, and scale 
radius, respectively. Given the cosmological paramet ers and haloes v irial m ass, we calculate 
these parameters following a method in Appendix of iNavarro et al.l ( 119971 ). 



The initial density profiles of the ICM are assumed to be those of a beta-model, 
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where Pgas,o and Tc are the central gas density and core radius, respectively. We assume that 
/ ? = 0.6, which is cons i stent with typical va lues obtained form the recent X-ray observations 
(IVikhlinin et al.ll2006l : ICroston et al.ll2008l ). On the other hand, it is not easy to obtain the 
observational relationship between Tc and r^ c onsidering diflScul t y of a ccurate determination 
of mass profile. Following the discussion in iRicker fc SarazinI ( J200ll ). therefore, we chose 
^c = ^s/2 as a fiducial value. The gas mass fraction is set to be 0.1 inside of the virial radius 
of each subcluster. It should be noted that the results will not be sensitive to the choice of 
this fraction unless the ICM dominates the DM in gravity. 

The velocity distribution of the DM particles is assumed to be an isotropic Maxwellian. 
The radial profiles of the DM velocity dispersion are calculated from the Jeans equation with 
spherical symmetry, so that the DM particles would be in virial equilibrium in the cluster 
potential of the DM and ICM, 



^(PDMCT^) 



GM{r) 



-PDMi 



(3) 



where a^ and M{r) are the one dimensional velocity dispersion at r and total mass inside r, 
respectively. G is the gravitational constant. We need boundary conditions to obtain a^(r) 



from differential equation ([3j), which is assume to be 



cr^ir^ir) 



where r^w is the virial radius of the halo. 
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The radial profiles of the ICM pressure are determined in a similar way so that the ICM 
would be in hydrostatic equilibrium within the cluster potential, 



dr 
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The boundary conditions for equation (JSD is 
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where ^^gpec is the specific energy ratio of the DM a nd gas at the virial radius. In case 
of the classical isothermal /?-model (see ISarazinlll986l ). /Jgpec should be equal to the fitting 
parameter (3 in equation (Ej) to obtain isothermal temperature profile. Although our model 
is diflFerent than the isothermal [3 model in underlying DM distribution, we assume /3spec=/? 
for simplicity. As a result, the temperature profile is not isothermal. 



2.3. Initial Conditions 



For typical runs, the DM masses within the virial radius of the larger and smaller 
subclusters are 5.0 x 10^^ Mq and 1.25 x 10^^ M©, respectively. Thus, the mass ratio is 4 : 1. 
It is useful to introduce an angular momentum parameter A = J|£'|-^/^/(GM^/^) in order to 
characterize oflF-center collisions, where J, £", and M are the angular momentum of the two 
subclusters around the center of masses, the binding energy between the two, and the total 
mass, respectively. Please note that A is the ratio between th e actual ang:ular velocity a nd 
the angular velocity needed to provide rotational support (see lBinnev fc Tremaind 119871 ). 



Given each subcluster's mass and an angular momentum parameter of the system, w e 
estimate "typical" initial conditions for cluster mergers in a similar way of lTakizawal (119991 ). 
How to estimate these conditions is described in detai l in Appendix, which is a natural 
generalization of the method for head-on collisions (see iTakizawal 1 19991 . section 2) to cases 
including off-center collisions. The coordinate system is taken in such a way that the center 
of mass is at rest at the origin. Two subclusters are initialized in the xy-plane, separated 
by a distance \/{Ri + R2Y — b'^ in the x-direction, and b in the ^-direction, where Ri^2 and 



b are the virial radii of the each subcluster and the impact parameter, respectively. The 
initial relative velocity is directed along the x-axis. The centers of the larger and smaller 
subclusters were initially located at the sides of x < and x > 0, respectively. 

At present, we have only limited information about magnetic field configurations in the 
intracluster space. Roughly speaking, however, they have random structures whose power 
spectrum can be approximated by a power-law in a smaller scale, and the mean field strength 
is an increasing function of the gas density. Taking these features into account, we make 
the initial magnetic field as follows. First, we generate random Gaussian vector potential 
in the wave- numbers space with a zero mean and single power-law spectrum, A{k) oc k~^. 
We assume /x = 5/3. Then, the generated potential is transformed into the real space using 
a three-dimensional Fast Fourier Transform and multiplied by a factor of pgas, which is 
expected assuming a uniform spherical collapse and fiux freezing. A random and divergence- 
free magnetic field is calculated via S = V x A. The normalization of the magnetic field 
is determined so that the total magnetic energy is one percent of the thermal one. Strictly 
speaking, each subcluster in the initial conditions is not in hydrostatic equilibrium with 
magnetic field. However, it does not a matter in our simulations because the magnetic 
energy is much less than the thermal one as we wrote above. 

The parameters for each model are summarized in tabled! where M^, r^^i^ Ci^ Ni are 
the total mass, scale radius, concentration parameter, and number of A^-body particles for 
the i-th subcluster, respectively. Agrd is the total number of grid points. RunST, which is a 
simple head-on collision case with the mass ratio of 4:1, is regard as a standard run. RunOC 
is a case of an off-center merger with A = 0.05. RunMM is a relatively minor merger case 
with the mass ratio of 8:1. RunLR is a lower resolution run to investigate how the results 
depend on numerical resolution. 



3. Results 

The left panels of Figure [1] show snapshots of the density (contours) and temperature 
(colors) on the x^-plain, which is perpendicular to the collision axis, at t = 1.11, 1.56, 2.00, 
2.89, and 3.78 Gyr of RunST. Right ones show the same but for the magnetic field strength. 
At t = 1.11 Gyr, two subclusters are approaching each other, and temperature between two 
density peaks is higher than elsewhere. However, magnetic field strength in this region is 
lower than that in the density peaks. This is simply because initial mean magnetic field 
strength increase with the density. Although the magnetic field there is slightly amplified 
by the adiabatic compression, this does not compensate the trend of the initial magnetic 
field. At t = 1.56 and 2.0 Gyr, slightly amplified magnetic field is seen just behind the 



shock. However, two kinds of structure appear more notably. One is a relatively strong 
magnetic field region associated with a contact discontinuity between the ICM originated 
from the smaller and larger subclusters, which will be recognized as a cold front with strong 
magnetic field in the observational point of view. As a result, a cool region surrounded by 
the magnetic field appears. The other is an ordered magnetic field in x-direction just behind 
the subcluster. This is because fiow behind subcluster is converging on the collision axis. As 
a result, the magnetic field are collected and amplified because of compression. At t = 2.89 
Gyr, the bow shock has gone outside of the panel, the field structure associated with the 
contact discontinuity and ordered field along the collision axis are still present. In addition, 
eddy-like field structures start to appear above and below the ordered field along the collision 
axis in x^-plain. These structures grow as Kelvin-Helmholtz instabilities develop and become 
clearer at t = 3.78 Gyr 

The initial conditions of RunST have an axial symmetric structure. Thus, the above- 
mentioned features could change significantly in off-center mergers. Figure [2] shows the same 
as figured! but for RunOC. Again, the prominent magnetic structures are associated with 
a contact discontinuity rather than a bow shock. Because of the asymmetry, the magnetic 
field at the contact discontinuity is stronger in the side closer to the lager cluster core (or 
lower side in each panel) at t = 2.00 and 2.89 Gyr. A cool region surrounded by the 
magnetic field certainly appears, but the field structure in the side farther to the larger 
cluster core (or upper side in each panel) becomes less clear at t = 3.78 Gyr. As in the 
case of RunST, an ordered and relatively strong magnetic field structure appears behind the 
moving substructure. However, this is not just behind the clump but shifted towards the 
lager cluster side. Eddy-like magnetic fields generated by the Kelvin-Helmholtz instabilities 
are clearer in the side further to the larger cluster. 

Cluster mergers cause various characteristic magnetic field configurations as we see 
before. Unfortunately, it is very difficult to observe the magnetic field structure itself directly 
at present. The observation of Faraday rotation measure is an indirect, but quite useful way 
to obta in information of the ma gnetic field structure. The rotation measure (RM) is given 



by (see lRybicki fc Lightmanlll979i ) 



RM[rad/m2] = 812 j n^[cm-^]B\\[iJiG\dl, (7) 

where n^ and B\\ are the electron number density and line-of-sight component of the magnetic 
field, respectively. Figure [3] presents snapshots of the RM map for RunST. Left and right 
panels show the maps seen from the z and x axis, respectively. Please remember that the 
collision axis is along the x-axis. When the line-of-sight is perpendicular to the collision 
axis, RM maps give us relatively rich information. Cool regions with high magnetic field are 
seen as those with high absolute values of RM. On the other hand, the structure is relatively 



featureless when we observe the system along the collision-axis. The absolute values of RM 
tend to be higher in the right panels, which is not surprising because the density distribution 
is elongated along the collision axis, and because there is the ordered magnetic field in the 
same direction behind the subclump. 

Figure m shows profiles along the collision axis at t = 2.00 Gyr of RunST for (a) pressure, 
(b) electron number density, (c) temperature, (d) x-component of velocity, (e) x-component 
of the magnetic field, and (d) absolute value of the magnetic field perpendicular to the axis. 
There is a bow shock at x c^ —1.5, where a clear discontinuity are seen in the profiles of 
pressure, density, temperature, and x-component of velocity. A somewhat blunt contact 
discontinuity is Sit x c^ —1^ where density and temperature have a jump, but pressure and 
x-velocity are smooth. The magnetic field perpendicular to the axis also increase across 
the bow shock. The tangential magnetic field is amplified between the shock and the con- 
tact discontinuity. As a result, the cool region is wrapped by the magnetized layers. This 
characteristic structure probably have an infiuence on the transport process such as heat 
conduction, which will be discussed later. 

Because minor mergers are more frequent in actual situations, it is interesting how the 
results change in minor mergers. Figure [5] shows the same as figure [4] but for RunMM, 
where mass ratio is 1:8. Mach number of the bow shock becomes lower, which is clearly seen 
in the profile of pressure, density, temperature, and x-component of velocity. As a result, 
amplification of the magnetic field component perpendicular to the collision axis around 
contact discontinuity is less prominent, where typical values are roughly halves of those in 
RunST. On the other hand, typical strength of the x-component of the magnetic field does 
not change very much. 

It is possible that our results, especially on the small scale structures around the cold 
front and bow shock, depend on numerical resolution. To assess this, we perform a lower 
resolution run as RunLR. Both the numbers of grid points and particles in RunLR are eighths 
of those in RunST. The initial conditions of RunLR is essentially the same as those of RunST 
though the largest wave- number components in the magnetic field are not included. Figure [6] 
shows the same as figured] but for RunLR. Although the overall features in pressure, density, 
temperature, and x-component of velocity are similar to those of RunST qualitatively, it is 
obvious that the width of the contact discontinuity is broader. Width of the bow shock is 
also slightly broader. As for magnetic field, although the small scale structures disappear, 
amplification of the field component perpendicular to the collision axis near the contact 
discontinuity is clearly seen again. 



4. Summary and Discussion 

We perform A^-body and MHD simulations of merging clusters of galaxies. We find that 
cluster mergers cause various characteristic magnetic field structures because of strong bulk 
fiow motion. The magnetic field component perpendicular to the collision axis is amplified 
especially between a bow shock and contact discontinuity. As a result, a cool region wrapped 
by the field lines appears. A relatively ordered field structure along the collision axis appears 
just behind the moving substructure. Eddy-like field structures are also generated by Kelvin- 
Helmholtz instabilities. Although the detailed structures change slightly similar features are 
seen in off-center mergers. RM maps have some information about the magnetic field struc- 
ture. The above-mentioned characteristic structures are partly recognized in the RM maps 
when we observe the merging systems in the direction perpendicular to the collision axis. 
On the other hand, RM maps observed nearly along the collision axis are less informative in 
this respect. Typical absolute values of the RM become higher when the system is observed 
along the collision axis because of both the density distribution elongated toward the axis 
and ordered magnetic field along the same direction. In minor mergers, amplification of the 
magnetic field component perpendicular to the collision axis becomes less prominent. The 
results about magnetic field along the collision axis does not change very much. Although 
numerical resolution effects have the significant impact on the small scale structures such as 
widths of a bow shock and contact discontinuity and detailed small scale fiuctuations of the 
magnetic field, our results about overall global structures are reasonably reliable. 

Some observational and theoretical studies suggest that heat conduction in the ICM is 
suppressed from the Spitzer value. For example, a sharp temperature jump across a cold 
front in A2142 requires that the heat conductivity is reduced by a factor o f between 250 



and 2500 at least in the direction across the front (JEttori fc FabianI l2000l ) . The spatial 



temperature variations in the central region of A754 also sugge st that the conductivit y 



is at least an order of magnitude lower than the Spitzer value ( JMarkevitch et al.l l2003l ) . 
It is well-known that fine-tuning of heat conductivity is necessary in order to reproduce 
observational natures of cool cores assuming that radiative c ooling in the cores balance s 



with the conductive heating from the outer part of the clusters (jZakamska fc Narayanll2003l ). 
Although the detailed process involved there is still unclear, magnetic fields likely play an 
important role. Heat conduction in the direction across the magnetic field lines is likely 
suppressed because electrons cannot travel freely in that direct i on. U sing two dimensional 



MHD simulations with anisotropic heat conduction, lAsai et al.l (120041 ) shows that a moving 
subclump naturally form the magnetic field structure along the contact discontinuity and 
that the temperature jump there is maintained if the heat conduction perpendicular to the 
field line is suffic iently suppressed. Similar results are obtained in three dimensional cases 



( JAsai et al.l 120071 ). Basically, our simulations confirmed their results about magnetic field 
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configurations in more realistic situations, though the heat conduction is not included. 

We show that Faraday rotation measure maps are useful to obtain information about 
characteristic magnetic field structures formed by cluster mergers. This could be another 
probe of internal dynamics of clusters. However, there is a serious problem in this method 
at present. We need polarized radio sources such as radio lobes and/or AGNs in and/or 
behind the clusters of galaxies to do that. In other words, we are able to obtain RM maps 
in the limited regions where we have suitable polarized sources by chance. Naturally, these 
do not always correspond to the regions that we are interested in. However, this difficulty 
could be overcame if we have suitable polarized sources tha t covers a cluster entirely. One 



possible solution is to use the CMB as the polarized source ( iQhno et al.ll2003l ). though this 
is still challenging in the present status of CMB observations. There are a lot of observations 
that are ongoing and planned for measuring the CMB polarization. We hope that future 
observations of the CMB polarizations will enable us to make RM maps that cover clusters 
entirely, which would give us an important clues to understand the internal dynamics as well 
as the magnetic field structures in clusters of galaxies. 

In actual clusters, the magnetic field have random components in small scales. It is very 
difficult to treat such scales by numerical simulations that concern the global structures. 
Clearly, observed RM maps are inffuenced by the small scale magnetic field fiuct nations, 
which are not considered in RM maps calculated from our results. As a result, coherent 
length of the magnetic fields tends to be overestimated effectively, which means absolute 
values of RM from the simulations are also overestimated. However, it is probable that 
global spatial patterns of RM maps are nearly independent of such small scale ffuctuations. 
In addition, it is also probably robust that higher RM are expected when the merging systems 
are observed along the collision axis. 
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M. T. was supported in part by a Grant-in- Aid from the Ministry of Education, Science, 
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A. An Estimation Method of the Initial Conditions of Mergers 

The estimation method introduced here is a natural extension of that for head-on col- 



hsions in iTakizawal (119991 ) to cases including off-center collisions. We consider a merger of 
two subclusters with masses Mi and M2 {Mi > M2). Let us consider a region that contains 
both subclusters. It expands following the cosmological expansion at first. However, the 
expansion is decelerated and then the region will collapse if the mean density is higher than 
the critical one. We assume that they are separated by the distance of 2rta at the maximum 
expansion epoch. In case of head-on collisions, it is obvious that they are at rest for their 
center of mass at that time. In case of off-center mergers, on the other hand, they have 
the relative tangential velocity '^0 at that time. Taking into account the conservation of 
energy and angular momentum between at the maximum expansion epoch and just before 
the merger, we obtain following equations, 

nu.M. 1 GM1M2 1,0 



2rta 2 ^ i?i + i?2 2 

2Mvon. = Mvb, (A2) 

where i?i, i?2, v^ and b are the virial radius for each subcluster, the initial collision velocity, 
and the initial impact parameter, respectively. M = MiM2/(Mi + M2) is the reduced mass 
of the system. It is natural that the virial radius has correlation with the mass. According 
to a spherical colla pse model, the virial mass is proportional to the square of the virial radius 



see 



Peebleslll980l) . Thus, we assume that the following relation. 



iJ,= (^) fi, (A3) 

The spherical collapse model also tell us that the turn around radius is twice of the virial 
radius. Thus, we also assume. 

Using equations (lAll) . flA2p . flASp . and (IA4I) . we eliminate Vq and obtain v as follows, 

(A6) 



where 

rv = 

Ml 



a - #. (A7) 
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It is useful to introduce an angular momentum parameter A = J\E\^^'^ /{GM^^'^) in order to 
characterize off-center mergers, where J, E, and M are the angular momentum of the two 
subclusters around the center of masses, the binding energy between the two, and the total 
mass, respectively. Please note that A is the ratio between th e actual angular velocity a nd 



the angular velocity needed to provide rotational support (see lBinnev fc Tremaind 119871 ). 
From equations flAip and flA2p . we obtain 



A = 



{GMiRif/^ (1 + a)7/2 1 1 + ai/3 ~ 2GM^{1 + a 



-.1/2 

} • (A8) 



Given Mi, i?i, a, and A, therefore, we are able to calculate v and h from equations flASp and 
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Fig. 1. — Left: Snapshots of the density (contours) and temperature (colors) distributions 
on the z = surface at t =1.11, 1.56, 2.00, 2.89, and 3.78 Gyr of RunST. Right: Same, but 
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Fig. 2. — Same as Fig. [U but for RunOC. 
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Fig. 3. — Snapshots of the Faraday rotation measure map of RunST at the same epochs as in 
Figdl Left and right panels show the rotation measure when the hne-of-sight is perpendicular 
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(b) Density 
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Fig. 4. — Profiles along the x-axis at t = 2.00 Gyr of RunST, for (a)pressure, (b)electron 
number density, (c)temperature, (d)x-component of velocity, (e)x-component of magnetic 
field, and (d)absolute value of magnetic field perpendicular to the x-axis, respectively. 
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Fig. 5. — Same as Fig. [H but for RunMM. 




(f) 




-2 -1.5 



-0.5 



[Mpc 



20 



(a) Pressure 



(b) Density 



(c) Temperature 




-1.5 -1 -0.5 




(e) 




-1.5 -1 -0.5 



X [Mpc] X [Mpc] 

Fig. 6. — Same as Fig. [H but for RunLR. 
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Table 1: Model Parameters 
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Parameters 


RunST 


RunOC 


RunMM 


RunLR 


Mi/M2(lOi4M0) 


5.0/1.25 


5.0/1.25 


5.0/0.625 


5.0/1.25 


^-8,1/^8,2 (kpc) 


256/139 


256/139 


256/104 


256/139 


C1/C2 


6.12/7.08 


6.12/7.08 


6.12/7.56 


6.12/7.08 


A 


0.0 


0.05 


0.0 


0.0 


iVl/iV2 


1677721/419431 


1677721/419431 


1864135/233017 


209715/52429 


A^grd 


(256)3 


(256)3 


(256)3 


(128)3 



